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We show numerically that in a binary system of Yukawa particles, a dispersity driven disordering 
transition occurs. In the presence of quenched disorder this disordering transition coincides with a 
marked increase in the depinning threshold, known as a peak effect. We find that the addition of 
poorly pinned particles can increase the overall pinning in the sample by increasing the amount of 
topological disorder present. If the quenched disorder is strong enough to create a significant amount 
of topological disorder in the monodisperse system, addition of a poorly pinned species generates 
further disorder but does not produce a peak in the depinning force. Our results indicate that for 
binary mixtures, optimal pinning occurs for topological defect fraction densities of 0.2 to 0.25. For 
defect densities below this range, the system retains orientational order. We determine the effect 
of the pinning density, strength, and radius on the depinning peak and find that the peak effect is 
more pronounced in weakly pinning systems. 

PACS numbers: 64.60.Cn,64.60.Ht,82.70.Dd,74.25.Qt 



I. INTRODUCTION 

In two dimensional (2D) systems, an ordered phase 
can undergo an amorphization transition in the presence 
of random quenched disorder. At zero temperature and 
for sufficiently strong quenched disorder, this transition is 
characterized by the appearance of isolated dislocations 1 . 
Disordering transitions can have a profound effect on 
transport properties. One of the best known examples 
of this phenomenon is the peak effect observed for vor- 
tices in type-II superconductors 2 -' 3 -' 4 * 5 -' 6 -. Here, an ordered 
vortex lattice is weakly pinned by random quenched dis- 
order; however, as the temperature or magnetic field is 
increased, the vortex lattice disorders and a sudden in- 
crease or peak in the pinning force occurs. The peak 
effect is known to be relevant to 2D and effectively 2D 
superconducting systems 3 ". 

An early explanation for the increase in the pinning 
force was that a disordered vortex system with numerous 
topological defects is much softer than an ordered vortex 
lattice, and as a result the vortices in the disordered state 
can shift positions in order to accommodate to the pin- 
ning landscape 2 * 6 -. A proliferation of topological defects 
has also been correlated with an increase in the effec- 
tive friction force in nanomaterials near the bulk melt- 
ing temperature where a peak effect type phenomenon 
can occur— Although the peak effect has been associ- 
ated with the appearance of topological defects, a com- 
prehensive understanding of exactly how the density or 
type of topological defects relates to the pinning effec- 
tiveness is still lacking. Interpreting the peak effect phe- 
nomenon is also complicated by thermal effects, and con- 
troversy remains over whether the disordering of the vor- 
tex lattice is predominately a melting phenomena 4 * 6 - or is 
nonthermal 9 . There are still several open issues related 
to the peak effect in superconductors and its connection 
to the shape of the current-voltage curves 1 -^, transient 
dynamics 1 * 1 -, and the magnetic field-temperature (H-T) 
phase diagram 12 . 



In two dimensional systems with binary or polydisperse 
interactions, a topological disordering transition occurs 
with increasing dispersit y 13 ' 14 : 15 . An open question is 
how quenched disorder affects a dispersity-driven amor- 
phization process and how the amorphization might alter 
the pinning effectiveness. A stiff lattice is poorly pinned 
by random disorder, and thus in a monodisperse assembly 
of repulsively interacting particles on quenched random 
disorder, the depinning threshold decreases when the lat- 
tice is stiffened by increasing the repulsion between the 
particles 1 -. A similar decrease of the depinning threshold 
with increasing particle-particle interactions occurs even 
when the monodisperse system contains some topologi- 
cal disorder and is no longer a perfect elastic lattice 16 . 
The situation may be different in a bidisperse system 
where the relative strength of the repulsive interactions 
between the two particle species can be independently 
tuned. If the interaction strength of one species is in- 
creased while that of the other species is held constant, 
the depinning threshold for motion over a random sub- 
strate may decrease; however, an amorphization of the 
bidisperse lattice occurs when the difference between the 
repulsive interactions of the two species is large enough. 
The effective softness of the lattice increases sharply once 
topological defects appear in the system, and thus the 
depinning threshold may increase rather than decrease 
when the particle interaction strength is increased across 
the amorphization transition. 

Polydisperse systems with quenched disorder can be 
used as models for exploring the interaction between de- 
fects and pinning since the number of topological de- 
fects can be controlled readily. Binary and polydis- 
perse particle interactions appear in a variety of sys- 
tems which can also have quenched disorder, including 
vortices in Bose-Einstein condensates**^, electron bub- 
ble mixtures*^, mixtures of Abrikosov and Josephson 
vortices 1 -^, and colloidal system o 20 ' 21 . In two dimensional 
colloidal systems, it was recently demonstrated that in- 
teractions with a very long screening length can be real- 
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ized experimentally^ and that disordered phases occur 
which are likely induced by polydispersitj*^. Two dimen- 
sional colloidal disordering transitions and pinning phe- 
nomena in the presence of quenched disorder have also 
been demonstrated experimentally^. Interaction disper- 
sity can also arise in the pinning of a nanofriction system. 
It would be very useful to understand how to control the 
defect density and pinning in these types of systems by 
manipulating the polydispersity in the particle interac- 
tions. 

In this paper we consider a two dimensional model of 
two species of particles interacting via a repulsive Yukawa 
potential both with and without quenched disorder. In 
the absence of quenched disorder, we find a well defined 
low temperature disordering transition as a function of 
dispersity. In the presence of weak quenched disorder, 
this disordering transition is accompanied by a sharp in- 
crease in the depinning force or a peak effect. The peak 
effect phenomenon we observe is nonthermal and occurs 
due to the proliferation of topological defects. Our re- 
sults indicate that there is a complex interplay between 
the density of topological defects and the depinning force. 
We identify where the peak effect phenomenon occurs as 
a function of pinning strength, density, radius, and inter- 
particle interaction strength. 



II. SIMULATION 



We simulate a two dimensional system of bidisperse 
particles in a sample of size L x = L y = 36 with peri- 
odic boundary conditions in the x and y directions. The 
particles interact via a Yukawa potential 
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where R^j) is the position of particle Rij = |Rj — 
R 3 |, Ci(j) is the charge of particle 1/k is the screen- 
ing length, and e is the dielectric constant. The inter- 
action force prefactor Z is set to Z = 1 unless other- 
wise noted. In this work we fix kappa at K = 2.0. We 
consider a binary arrangement of Na and Nb particles 
of species A and B, respectively, which have different 
charges denoted by Ca and C'b- For a monodisperse 
system Cb/Ca = 1-0. The total number of particles 
N = Na + Nb and the density n = N/L x L y . In gen- 
eral we fix N — 856 and Ca — 1 and vary either the 
ratio Nb/N or Cb/Ca- The particles also interact with 
a random quenched background which is modeled as N p 
randomly placed parabolic traps of radius r p — 0.3, den- 
sity n p = N p /L x Ly = 0.66 and maximum force f p . The 
parabolic trap potential is similar to that used in previ- 
ous work2£, Vp(Rik) = -(f p /2r p )(Rik - r p f for R ik < r p 
and zero interaction for > r p , where Rik — |R» — Rjt| 
and Rfe is the position of trap k. The particles evolve 
under Brownian dynamics, performed by integrating the 
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FIG. 1: (a) Fraction of defected particles Pd vs the ratio of 
the particle interaction strength Cb/Ca for a mixture with 
Nb/N = 0.5, Ca = 1.0, and f p = at different net particle 
densities n = 0.97 (filled triangles), n = 0.66 (filled circles), 
and n = 0.32 (open squares), (b) Phase diagram of n vs 
Cb/Ca for the same system. 



overdamped equation of motion 
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where r\ is the damping constant. The external dc driv- 
ing force Fd = F&X is slowly increased from zero in 
increments of<5Fc = 2 x 10~ 5 applied every 2.5 x 10 4 
simulation steps. We have found that slower increment 
rates do not change the results. The thermal force F T 
is modeled as random Langevin kicks with (F T ) = 
and (F T (t)F T (£')) = 2rjk B T5{t - t'). Since we are inter- 
ested in nonthermal effects, we consider a low tempera- 
ture T/T m = 0.15, where T m is the melting temperature 
of a monodisperse system with Nb/N = at density 
n = 0.66. The initial particle configurations are obtained 
using two techniques which produce identical results. In 
the first, we perform simulated annealing, while in the 
second, we place the particles in a triangular lattice, and 
when the dispersity is large enough defects naturally ap- 
pear. 



III. EFFECT OF PARTICLE DISPERSITY 
WITHOUT QUENCHED DISORDER 

To demonstrate that this system exhibits a dispersity 
driven disordering transition in the absence of quenched 
disorder, we consider a sample with Nb /N — 0.5 and 
Ca — 1.0 at different particle densities n = 0.97, 0.66, 
and 0.32 for f p = 0. In Fig. 1(a) we plot the fraction 

of defected particles Pd — N^ 1 5^ i=1 [l — 5(6 — z,)] as a 
function of Cb /Ca at the three different densities, where 
Zi is the coordination number of particle i obtained from 
a Delaunay triangulation. In each case, when Cb/Ca is 
near 1, the system forms a triangular lattice free of topo- 
logical defects and Pd = 0. Once the dispersity is strong 
enough, a disordering transition occurs with a prolifer- 
ation of defects and Pd rises above zero. The dominant 
type of defects we observe are fivefold and sevenfold coor- 
dinated particles, although in the very disordered states 
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FIG. 2: (a,b) Particle positions in a 24 x 24 portion of a 
sample with Nb/N — 0.5, n = 0.66, and f p = 0. Large 
circles: species A; small circles: species B. (a) The or- 
dered regime at Cb/Ca = 1-2. (b) The disordered regime 
at Cb/Ca = 2.0. (c,d) The corresponding Delaunay trian- 
gulations at (c) Cb/Ca = 1.2 and (d) C s /Ca = 2.0. Dark 
circles: fivefold coordinated particles; light circles: sevenfold 
coordinated particles. 



it is possible to find a small fraction of fourfold or eight- 
fold coordinated particles. As the density of the system 
increases, the dispersity Cb/Ca at which the disorder- 
ing transition occurs decreases since the particle-particle 
interactions are stronger at higher density and thus the 
strain energy induced by the dispersity in the particle 
interactions increases. 

By conducting a series of simulations we map the 
order-disorder transition line as a function of density n 
and polydispersity Cb/Ca, shown in Fig. 1(b). The or- 
der to disorder transition is defined as occurring when 
the fraction of defected particles reaches Pd = 0.2. If a 
different cutoff value of Pd is used to identify the tran- 
sition, the general features of the phase diagram are un- 
changed but the precise location of the transition can 
shift. At low densities n < 0.06 the system is in a liquid 
state even for Cb/Ca = 1.0 since the particle-particle 
interactions are not strong enough to overcome the ther- 
mal forces. Since the simulation is performed at finite 
temperature, the thermal forces become important when 
the particle-particle interactions are weak enough. For 
a system with a finite particle-particle interaction range, 
the particle density can be decreased to a point at which 
there is no elasticity since adjacent particles are no longer 



interacting on average. This is what gives rise to the low 
density liquid state. 

In Fig. 2(a) we illustrate the particle positions for 
a system at n = 0.66 with Cb/Ca = 1.2 in the or- 
dered regime. The corresponding Delaunay triangulation 
shown in Fig. 2(c) indicates that all of the particles are 
sixfold coordinated and form a triangular lattice. Figure 
2(b) shows the particle positions in the same system at 
Cb/Ca — 2.0 in the disordered regime, while Fig. 2(d) 
illustrates that in this regime the lattice is filled with 
fivefold and sevenfold coordinated defects. 

In the phase diagram of Fig. 1(b), we only distinguish 
between the topologically ordered and topologically dis- 
ordered states, and indicate the reentrant disordering 
transition into a liquid state that occurs at low den- 
sities. At intermediate densities, the disordered phase 
may have glasslike features, implying that there could be 
another line on the phase diagram between a liquidlike 
and glasslike state. Determining whether such a line is 
present or absent is beyond the scope of the present pa- 
per. It is likely that the actual lowest energy state for 
the polydisperse system is phase separated; however, we 
have never observed such a state. Recent experiments 
and simulations of binary colloidal systems with disper- 
sity similar to what we consider here also produced no 
phase separated states, but did show some evidence for 
clustering 20 . In simulations by Sadr-Lahijany et al^, an 
intermediate hexatic phase appeared in certain regions of 
the density-polydispersity phase diagram. The hexatic 
phase is characterized by an algebraic decay in the orien- 
tational correlation function. In general, hexatic phases 
are difficult to observe and there are still open questions 
about the nature of this phase. It is beyond the scope of 
this work to address the possible existence of a hexatic 
phase in clean systems; instead, we focus on the regimes 
with quenched disorder. 



IV. EFFECT OF QUENCHED DISORDER 

In Fig. 3(b) we plot Pd vs Cb/Ca for a mixture with 
Nb/N — 0.5 and n = 0.66 in the presence of a ran- 
dom pinning potential with f p = 0.1 and n p = 0.66. At 
this pinning strength the system is free of dislocations 
for Cb/Ca < 1-6. Near Cb/Ca — 1-75 an order-disorder 
transition occurs. In Fig. 3(a) we plot the correspond- 
ing critical depinning force F c which is the value of Fjj 
at which the velocity of the particles in the direction of 
the drive exceeds 5 x 10 -4 . In general there are two 
predominant pinning regimes. For the defect-free lattice 
at Cb/Ca < 1.6, the depinning occurs elastically with- 
out the further generation of defects. In this case the 
depinning threshold is well defined as the entire lattice 
moves at the same average velocity. Once defects begin 
to appear around Cb/Ca = 1.75, the depinning becomes 
plastic and a portion of the particles can be moving while 
another portion remains pinned, leading to very inho- 
mogeneous and intermittent velocity bursts. The large 
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FIG. 3: (a) Depinning force F c and (b) fraction of de- 
fected particles Pd vs Cb/Ca for a mixture at n = 0.66 with 
Nb/N = 0.5 on a random pinning substrate with f p — 0.1 
and n p = 0.66. (c,d) F c and Pd versus Cb/Ca for a system 
with the same parameters except with f p = 0.4. 



fluctuations associated with plastic depinning can also 
lead to long transient responses where a portion of the 
system is initially moving when a drive is applied, but 
over time the entire system becomes pinned. In order 
to accurately measure a depinning threshold in the topo- 
logically disordered regime, sufficiently long run times 
at each force increment are needed in order to obtain a 
smoother velocity- force curve. We also note that chang- 
ing the value of the threshold velocity for the depinning 
measurement does not change the general features of the 
results. 

Figure 3(a) shows that the depinning force F c initially 
decreases with increasing Cb/Ca in the ordered regime, 
but then starts to shift upward once Cb/Ca = 1.2. At 
the order-disorder transition which occurs at Cb/Ca = 
1.75 there is a sharp increase in F c followed by a slow de- 
crease for increasing Cb/Ca- We term this phenomenon 
a dispersity driven peak effect. This result demonstrates 
that a peak effect can arise via a completely nonthermal 
disordering process. 

For a monodisperse system with Cb/Ca = 1, increas- 
ing the magnitude of the repulsive interaction Cb always 
reduces i 7 !—. In contrast, the results in Fig. 3(a) show 
that increasing the charge of only a fraction of the parti- 
cles can actually increase the depinning force. The max- 
imum in F c occurs when a fraction Pd = 0.2 of the par- 
ticles are defected. As Cb/Ca increases further above 
Cb/Ca — 1.75, Pd increases; however, the depinning 
force decreases. This result indicates that although a 
peak effect occurs at the onset of defect proliferation, 
simply having more topological defects does not directly 
translate to a higher depinning force. 

For a monodisperse system, increasing the strength 
of the repulsive particle-particle interactions decreases 
the depinning threshold. In the bidisperse system in 



Fig. 3(a), when Cb/Ca is increased above 1 by only a 
small amount, the net distortion of the lattice is small 
and no dislocations are induced since the system is still 
nearly monodisperse. Thus in this regime, increasing 
Cb/Ca increases the overall effective lattice stiffness and 
causes the depinning threshold to decrease. Once Cb / Ca 
is large enough [Cb/Ca > 1.2 in Fig. 3(a)], more signifi- 
cant distortions of the lattice occur due to the fact that 
the two particle species are randomly interspersed, and 
these distortions effectively soften the lattice and increase 
the depinning threshold. As Cb/Ca continues to in- 
crease, the distortions become large enough to induce the 
formation of dislocations [Cb/Ca = 1-75 in Fig. 3(a)], 
and the order-disorder transition occurs. The depinning 
threshold increases rapidly since the particles can shift 
into optimal pinning locations once the dislocations ap- 
pear, and the maximum in F c occurs at Cb/Ca = 1-8. 
When Cb/Ca is increased above 1.8, defects continue to 
proliferate in the lattice but the particles arc already in 
their optimal pinning locations so no further enhance- 
ment of the depinning threshold occurs. Instead, as the 
particle-particle interaction force becomes stronger, the 
local stiffness of the lattice increases, shifting the parti- 
cles away from the optimal pinning locations and causing 
the depinning threshold to decrease again with increasing 
Cb/Ca > 1.8. 

To show that the order-disorder transition is respon- 
sible for the peak in the depinning force, in Fig. 3(c,d) 
we plot F c and Pd vs Cb/Ca for a system with the same 
parameters as in Fig. 3(a,b) but with a stronger pinning 
force of f p = 0.4. For this value of f p , the monodisperse 
system Cb/Ca = 1.0 already contains a significant frac- 
tion of topological defects, with Pd = 0.34. As Cb/Ca 
increases from 1 the defect density Pd decreases; how- 
ever, for Cb/Ca > 1.5 Pd begins to increase again. In 
contrast, the depinning force F c decreases monotonically 
for all Cb/Ca > 1.0. We have performed a series of sim- 
ulations for other values of f p and pinning densities n p 
and find the following general features. (1) A peak effect 
phenomenon occurs whenever there is an order to disor- 
der transition. (2) The peak in the depinning force occurs 
when a fraction of about Pd — 0.15 to 0.3 of the particles 
are defected. (3) If the system is already strongly disor- 
dered for the monodisperse case Cb / Ca = 1 , there is no 
enhancement in depinning force with increasing Cb/ Ca 
even when the fraction of topological defects increases. 



V. CHANGING THE PARTICLE SPECIES 
RATIO 

We next consider the effect of varying the species ra- 
tio Nb/N at fixed particle density n — 0.66 to ex- 
plore the disordering and depinning effects as the mix- 
ture varies from monodisperse species A, Nb/N — 0, to 
monodisperse species P>, Nb/N = 1. In Fig. 4(a) we 
plot the depinning force F c vs Nb /N for a system with 
f p = 0.1, n p = 0.66, and Ca = 1 at Cb/Ca = 1.5 
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FIG. 4: (a) The depinning force F c for varied species ratio 
N B /N at fixed n = 0.66, n p = 0.66, and f p = 0.1. Filled 
squares: Cb/Ca = 1.5. Open circles: Cb/C'a = 1.825. (b) 
The corresponding fraction of defected particles Pd- (c) F c 
vs Nb/N for the same set of parameters in (a) but with 
Cb/Ca = 2.5. Open circles: f p = 0.1. Filled squares: 
f p = 0.4 (this curve was divided by 10 for presentation pur- 
poses), (d) The corresponding P d curves. 



and Cb/Ca = 1.825. Figure 4(b) shows the corre- 
sponding P d curves. For Cb/Ca = 1.5, there is no dis- 
location induced transition at any fraction Nb/N and 
F c shows no enhancement but merely decreases above 
N B /N = 0.5. The depinning force F c = 0.003 for the 
monodisperse B system at Nb/N = 1 is less than that 
of the pure A system, F c = 0.005 at N B /N = 0, as 
expected for monodisperse systems due to the increased 
repulsive force between species B particles compared to 
species A. The depinning force for the pure B system is 
lower at Cb/Ca — 1-825 than at Cb/Ca — 1-5; however, 
for intermediate values of Nb/N in the Cb/Ca — 1-825 
sample there is a strong enhancement of F c over the de- 
pinning force for either the pure A or pure B systems. 
The enhancement of F c for Cb/Ca = 1.825 is associated 
with the creation of defects, as shown by the behavior of 
Pd in Fig. 4(b). The defect density reaches a maximum 
value of P d = 0.21 at N B /N = 0.5. The peak value of 
F c also falls at Nb/N = 0.5. These results show that, 
due to induced disorder, the effective pinning force for 
mixed species can be higher than that of cither of the 
pure species. 

In Fig. 4(c) we plot F c vs Nb/N for a system with 
Cb/Ca = 2.5 and C A = 1.0 at f p = 0.1 and f p = 0.4. 
For f p = 0.1, the system becomes increasingly disordered 
as Nb /N increases from zero until Pd reaches a maximum 
of P d = 0.48 at N B /N = 0.44, as shown in Fig. 4(d). 
For N B /N > 0.44, P d decreases back to P d = for the 
pure B system. There is a peak in F c at Nb/N = 0.2, 
where P d = 0.3. In general we find that for increasing 
Cb I Ca and fixed C A , the peak in F c shifts to lower val- 
ues of Nb/N. The peak in F c occurs when a fraction 
P d ~ 0.2 of the particles are disordered, and since higher 
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FIG. 5: (a) The depinning force F c for varied species ratio 
Nb/N at fixed n = 0.66, n p = 0.66, f p = 0.1, and Cb/Ca = 
1.5. Open diamonds: T/T m = 0. Filled circles: T/T m = 0.15. 



values of Cb /Ca disorder the system more effectively, a 
lower fraction of species B is needed to reach P d = 0.2 
as Cb/Ca increases. This result suggests that strong en- 
hancement of the pinning can be achieved by adding a 
few strongly repulsive particles to a pure system. For the 
case of f p = 0.4 in Fig. 4(c,d), the pure A case is already 
strongly defected, so as Nb/N is increased from zero, F c 
monotonically decreases while P d increases slightly and 
then decreases. 

We find that for decreasing f p , the maximum value 
pmax Q f j n |_ ne m j xec j sample increases relative to the 
value FP ure of F c in the pure samples with N B /N = 
or 1. This can be understood by a simple argument. For 
the pure samples, the pinning is collective and F^ ure cx 
F p , while at the mixed sample peak the pinning behaves 
more like single particle pinning with F™ ax « F p , so that 
F™ ax /FP ure oc l/F p . This indicates that in pure systems 
with weak pinning, the addition of a second species could 
have a very significant effect on the pinning properties. 



The results up to this point were obtained at a finite 
but low temperature. Most experiments with colloids 
are performed in regimes where some Brownian motion 
occurs, and our results show that the peak effect phe- 
nomenon is robust against the addition of a small tem- 
perature and thus could be observed in colloidal experi- 
ments. At higher temperatures, a significant amount of 
creep occurs in the disordered regimes and makes ob- 
taining an accurate depinning threshold difficult. This 
type of thermally activated motion will be described else- 
where. To illustrate that our results remain unchanged 
for T = 0, in Fig. 5 we plot F c vs Nb/N for the same 
system in Fig. 4(a) at Cb/Ca — 1-5, showing that there 
is a negligible difference between the T = and finite T 
results. 
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FIG. 6: (a) The depinning force F c vs pinning density n p for 
a system with fixed f p — 0.1, r p = 0.3, and Nb/N = 0.5 for 
the bidisperse system Cb/Ca = 1-85 (black circles) and the 
monodisperse system Cb/Ca = 1.0 (open squares), (b) The 
ratio R of F c for the bidisperse and monodisperse samples 
in (a), (c) The depinning force F c vs pinning radius r p for 
a fixed pinning density of n p = 0.66 at Nb/N = 0.5 for 
Cb/Ca = 1-85 (black circles) and Cb/Ca = 1.0 (squares), 
(d) The ratio R of F c for the bidisperse and monodisperse 
samples in (c). 



FIG. 7: (a) The depinning force F c vs pinning strength f p 
for a system with fixed n p — 0.66, n = 0.66, r p = 0.3, and 
N B /N = 0.5. Filled circles: Cb/Ca = 1.825; open squares: 
Cb/Ca = 1-0. (b) The ratio R of the bidisperse and monodis- 
perse critical depinning forces from (a), (c) The depinning 
force F c vs particle interaction strength Z for a system with 
fixed n p = 0.66, n = 0.66, r p = 0.3, and N B /N = 0.5. Filled 
circles: Cb/Ca = 1.825; open squares: Cb/Ca = 1-0. (d) 
The ratio R of the critical depinning forces from (c) . 



VI. VARIED PINNING DENSITY AND 
PINNING RADIUS 

We next vary several parameters of the pinning sites 
in order to understand how general the peak effect phe- 
nomenon is in binary systems with quenched disorder. 
In Fig. 6(a) we show F c vs pinning density n p for a 
sample with Nb/N = 0.5, n = 0.66, n p = 0.66, and 
f p = 0.1 at Cb/Ca — 1-825, corresponding to the 
peak value of F c in Fig. 3(a), and Cb/Ca = 1-0, cor- 
responding to a monodisperse system. To indicate the 
magnitude of the peak effect, we calculate the ratio of 
the bidisperse and monodisperse depinning thresholds, 
R = F C {C B /C A = 1.825) /F C {C B /C A = 1), and plot the 
result in Fig. 6(b). Although F c increases monotonically 
with rip for both values of Cb/Ca, F c is always higher 
in the polydisperse system. The ratio R goes through a 
maximum near n p — 0.66 when the density of particles 
equals the density of pins. At low n p , the bidispersity 
leads to little enhancement of the depinning threshold 
compared to the monodisperse system. When there are 
few pinning sites, most pinning sites can be occupied even 
in the monodisperse case, so softening the lattice by mak- 
ing the particles bidisperse and introducing dislocations 
does not create a significantly larger number of pinned 
particles. As a result, the enhancement of F c by the 
polydispersity is weak in this regime. Similarly, at high 
pinning densities there are so many pins available that 
less distortion of the lattice is required to permit most of 



the particles to occupy pinning sites, so the introduction 
of dislocations due to bidispersity does not create as large 
of an enhancement in F c compared to the intermediate 
pinning density regime. 

In Fig. 6(c) we plot F c versus pinning radius r p for 
a system with fixed n p = 0.66, n = 0.66, Nb/N = 0.5 
and f p = 0.1 at C B /C A = 1.825 and C B /C A = 1.0. We 
plot the corresponding ratio R of the depinning forces 
in Fig. 6(d). Here, R decreases monotonically with in- 
creasing r v . As the size of the pinning sites increases, 
the amount of lattice distortion required to permit most 
particles to occupy a pinning site decreases since a larger 
area of the sample is covered by the pinning sites. Thus, 
there is a decreasing enhancement of F c due to bidisper- 
sity as r p increases. 



VII. VARIED PINNING STRENGTH AND 
PARTICLE INTERACTION STRENGTH 

In Fig. 7(a) we plot F c vs pinning strength f p for a sys- 
tem with n p — 0.66, n = 0.66, r p — 0.3, and Nb/N = 0.5 
for the bidisperse case Cb/Ca = 1.825 and the monodis- 
perse case Cb/Ca = 1.0. The corresponding ratio R of 
the bidisperse to the monodisperse depinning threshold is 
plotted in Fig. 7(b). For large f p , both the monodisperse 
and the bidisperse systems are considerably defected and 
flwl, In this regime the critical depinning force F c de- 
pends linearly on f p , as expected for single particle pin- 
ning behavior. For f p < 0.2, the monodisperse system 



is defect-free and the depinning crosses over to collective 
pinning behavior with F c oc f p . In the same regime, true 
collective behavior does not occur for the bidisperse sys- 
tem and the depinning force falls off less steeply than f£, 
as seen in Fig. 7(a). Thus the ratio R between the bidis- 
perse and monodisperse depinning thresholds grows with 
decreasing f p , as shown in Fig. 7(b) for f p < 0.2. This 
result indicates that for weaker pinning, the peak effect 
phenomenon is more pronounced. 

To vary the overall particle-particle interaction 
strength, we adjust the value of the interaction prefac- 
tor Z in Eq. (1). For larger Z, the particles are more 
repulsive and the lattice is stiffer. We plot the effect 
of changing Z on F c in Fig. 7(c) for both monodisperse 
(Cb/Ca = 1) and bidisperse (Cb/Ca = 1.825) sam- 
ples with rip — 0.66, n = 0.66, r p — 0.3, f p — 0.1, and 
Nb/N = 0.5. The corresponding ratio R of the depin- 
ning threshold for the bidisperse and monodisperse sam- 
ples is shown in Fig. 7(d). We note that changing Z has 
an effect similar to changing k since the lattice becomes 
stiffer for smaller n. For Z > 1, the monodisperse system 
at C'b jC a — 1 is defect free, the depinning is elastic, and 
R increases with increasing Z . In this case, the effect of 
increasing Z is similar to the effect of decreasing f p . For 
Z < 1, the monodisperse system becomes topologically 
disordered and R approaches 1.0 as Z decreases. 

These results show that the peak effect phenomenon is 
observable in a polydisperse system whenever the effec- 
tive quenched disorder is weak enough that the monodis- 
perse system is free of defects. The monodisperse lattice 
is pinned collectively by the random disorder, and the col- 
lective pinning is always weaker than the single particle 
pinning which occurs in a topologically disordered poly- 
disperse particle lattice. In regimes where the quenched 
disorder is strong enough that both the monodisperse and 
the polydisperse systems are topologically disordered, 
there is little difference in the depinning force between the 
polydisperse and the monodisperse cases, and no peak in 
the depinning force appears. 



VIII. DEFECT DENSITY AND 
ORIENTATIONAL ORDER 

A general trend visible in Fig. 3(a,b) and Fig. 4(a,b,c,d) 
is that the peak in F c corresponds to a defect density of 
around Pd « 0.2. For defect densities Pd > 0.2, the de- 
pinning force F c decreases with increasing Pd- Recent 
simulations of two-dimensional monodisperse Yukawa 
particle systems have shown that as a function of defect 
density, the orientational correlations are lost when the 
fraction of defected particles is Pd > 0.22&. The orien- 
tational correlations are measured by means of the local 
bond orientational order parameter—, 

^(R,) = ^Ee 6 ^y (3) 
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FIG. 8: The value of the bond orientational correlation func- 
tion ge(r) at r = 5.0 versus the fraction of defected particles 
Pd in a system with Cb/Ca = 2.5 and f p = 0.1 at various 
particle ratios Nb/N. 



Here Ni is the number of nearest neighbors of parti- 
cle i and Oij is the angle between an arbitrary fixed 
reference axis and the bond connecting particles i and 
j. The bond orientational correlation function ge(r) = 
(iPq(R/)iPq(H' — r))/g(r), where the pair distribution 
function g{r) = (<5(R')<5(R' — r)). Orientational order is 
present when g§(r) decays algebraically with r. 

We measure ge(r) as a function of defect density Pd in 
a system with Cb/Ca = 2.5 and f p — 0.1 over a range 
of defect fractions from Pd — 0.0 to Pd — 0.5 obtained 
by varying Nb/N, as indicated in Fig. 4(d). In Fig. 8 
we plot the value of g§(r) at r = 5.0 versus Pd- The ori- 
entational order decreases monotonically with increasing 
Pd- We also find that ge (r) decays exponentially with r 
for Pd > 0.2, indicating that the system has lost orienta- 
tional order above Pd — 0.2. 

For defect densities below Pd = 0.2, the system is ori- 
entationally ordered and the depinning is elastic in na- 
ture. The defects soften the lattice locally, allowing the 
particles to become better pinned. As more defects are 
added and Pd increases, the lattice continues to soften 
until the orientational order is lost at Pd — 0.2, which 
corresponds to the elastic-plastic depinning transition. 
Once in the plastic depinning regime, the addition of 
more particles or more defects does not further increase 
the softness of the lattice. In our system, species B is 
more highly charged than species A and thus species B 
is in general less well pinned due to the higher strength of 
the particle-particle interaction compared to the particle- 
pin interaction. As Nb/N increases from zero, the de- 
fect density of the entire system increases, which would 
normally lead to more effective pinning, but clumps of 
species B particles can form and create localized poorly 
pinned areas, which may decrease the overall pinning ef- 
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fectiveness. Although it is beyond the scope of this pa- 
per to address, an open question is to understand ex- 
actly why a fraction of Pd = 0.2 defects corresponds to 
the loss of orientational order in two-dimensional sys- 
tems and whether this fraction Pd = 0.2 of defects is 
universal to all two-dimensional systems at the orienta- 
tional disordering transition regardless of the form of the 
particle-particle interactions. 

IX. SUMMARY 

To summarize, we have shown that for a binary as- 
sembly of repulsively interacting Yukawa particles, there 
can be an order to disorder dispersity driven transition 
induced by the addition of a more strongly repulsive 
species. In the presence of weak quenched disorder, the 
dispersity driven transition coincides with a peak effect 



in the depinning force. The depinning force for a mixture 
of particles can be higher than that for either of the pure 
species. If the quenched disorder is strong and the pure 
species sample is already disordered, the addition of a 
more strongly repulsive species monotonically decreases 
the depinning force. For a completely ordered system, if 
the addition of a new species does not induce dislocations, 
there is no enhancement of the depinning force. We also 
find that the optimal pinning occurs when the fraction 
of non-sixfold coordinated particles is around 0.2, which 
correlates with the loss of orientational order. Our results 
suggest that a peak effect phenomenon can occur with- 
out thermal fluctuations and that the effectiveness of the 
pinning can be tuned by adjusting the particle mixture. 

This work was carried out under the auspices of the 
NNSA of the U.S. DoE at LANL under Contract No. 
DE-AC52-06NA25396. 
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